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AND BOLTED-JOINT PROBLEMS 
By Robert D. Rosa 


SUMMARY 


The analogy "between the distritmtion of stresses in flat 
stiffened panels and the distribution of electric current in a 
ladder-type resistance network is used as the theoretical "basis 
of an electrical computer for the rapid solution of shear-lag 
problems. The computer, consisting of variable resistors and 
multiple -current sources, is described; and typical examples are 
given of its use. The analogy is extended to include bolted- 
joint problems, and an example is given also. 


INTRODUCTION 


When a load is applied along the flange at tne edge of a 
flat stiffened panel (fig. 1), some of the load is transmitted 
to the various stiffeners through the medium of the sheet, which 
is placed in shear by the action. Similarly, when two plates 
are connected by a bolted strap (fig. 2), the load is transmitted 
from the plates to the strap through the medium of the bolts. As 
shown in reference 1 and in appendix A, an analogy can be drawn 
between the distribution of forces in the cases previously 
described and the distribution of electric current in a ladder 
network. This analogy is very useful in obtaining a numerical 
solution, particularly in complicated cases which are not readily 
amenable to mathematical solution. The time saved in such cases 
is of appreciable importance. 

In reference 1 the solution of a shear-lag problem was 
obtained by means of an analogous network consisting of properly 
chosen fixed resistors. The power supply consisted of a single 
battery; the proper input currents were obtained from a voltage 
divider by trial and error. 

The present paper describes a pilot model of a more elaborate 
setup. The model uses variable resistors and several sources of 
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current, "both of which can "be introduced independently into the 
network. Being readily adaptable to the solution of more than 
one problem, this construction converts the apparatus into a 
computer . 

SYMBOLS 

A cross-sectional area, square inches 

B load carried by an individual bolt, pounds 
C bolt constant, inches per kip 

D diameter of bolt, inches 

E Young’s modulus, lcsi 

G shear modulus, ksi 

I longitudinal current, milliamperes 

K axial flexibility, inches per kip 

L axial length, inches 

P applied load, pounds 

P longitudinal resistance, ohms 

S shear force, po’uids 

V electric potential, millivolts 

b transverse width, inches 

h depth of beam, inches 

i current in transverse resistors, milliamperes 

p pitch, inches 

r resistance of transve - 3© resistors, ohms 

t thickness, inches 
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u axial displacement, inches 

a proportionality factor for resistance 

P proportionality factor for load 

a direct (normal) stress, psi 

t shear stress, psi 

Subscripts: 
b bolt 

F flange 

L longitudinal (stringer) 

m a particular transverse row of resistors or bolts 

n any transverse row of resistors or bolts 

p plate 

s sheet or strap 

T total 

av average 

DESCRIPTION OF COMPUTER 


The front panel of the computer supports the variable 
resistors constituting the network, together with the resistor 
control knobs and dials, as shown in figure 3* In order to set 
any element to a given resistance, the dial is turned to the 
proper reading as determined by the calibration constant for 
that element since the dials are graduated in 100 arbitrary 
divisions rather than in ohms. The dial constants were individually 
determined; consequently, it was possible to set the resistance 
elements to tl percent or i2 ohms, whichever unit was the greater- 

The resistor elements are connected to form a ladder-type 
network as shown in the wiring diagram in figure 4. Each element 
has a maximum and minimum resistance of approximately 1000 ohms 
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and 10 ohms, respectively. Turning the dial to the extreme low 
end short-circuits the minimum resistance of the element; whereas 
turning the dial to the high end open-circuits the element after 
the 1000 -ohm limit is reached. 

Current may he independently introduced into the network at 
eight points along the two sides opposite the ground connections. 
The current is supplied hy electronic regulating circuits (fig. l), 
which are designed to give negligible reaction upon each ether 
and to he independent of the resistance setting of the network. 

The effect of line-voltage variation is also negligible . The 
output of the regulating circuits can he varied from about 1 
to il milliamperes hy means of control dials on the front of the 
panel . Figure 6 is a rear view of the computer and shows the 
arrangement of the resistors and the current supplies. 

The current in any element, as well as the output of the 
electronic supply units, is determined hy plugging in a meter 
to a current -measuring jack associated with the desired element. 

A high-grade mi 11 i ammeter is used and allows the current to he 
measured with an accuracy of about l/U percent of full scale. 

The chief sources of error in the apparatus itself arise from 
inaccuracy in resistor values and in the current measurement. In 
the present model, an attempt was made to make a reasonable com- 
promise between absolute minimization of all sources of error on 
one hand and cost and convenience on the other hand. For instance, 
the variable resistors are of comparatively high resistance; thus, 
the errors due to variable contact resistance and meter- insertion 
loss are reduced. If the apparatus is allowed to warm up for 
about 11 minutes, errors due to further heating are negligible. 

The errors inherent in the method, when used for shear-lag 
problems, are d_ue to the division of the structure into a finite 
number of bays (corresponding to an electrical transmission line 
with lumped constants) and due to the inclusion of the sheet area 
with the stringer a’-ea to account for the axial stress in the 
sheet. These errors have been discussed and evaluated in 
reference 1. 


ANALOGY FOR SHEAR -LAG PROBLEMS 


As shown in reference 1 the correspondence between the 
mechanical structure and the electrical network is given in the 
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notation of the present paper "by the following table: 


Structure 

Network 

Quant i ty 

Symbol 

Quantity 

Symbol 

Axial load 

P 

Longitudinal current 

I 

Axial flexibility 

w 

il 

gif 

Longitudinal resistance 

1 

B 

Shear force 

S 

J Transverse current 

i 

Shear flexibility 

b 

GtL 

Transverse resistance 

r 

Axial displacement 

u 

Potential 

V 


Vritten in equation form with arbitrary constants 3 and a, 
the following relations hold: 


I = 3P = pcrA 
i * 3S = 3rLt 



( 1 ) 

( 2 ) 

(3) 

(M 

(5) 


The boundary conditions are not included in these equations and, 
therefore, have to be determined from the individual problem, in 
which current corresponds to force and potential differences correspond 
to displacements. At a fixed boundary no displacement occurs; 
hence, there must be no potential difference along that part of the 
network representing the boundary. The electrical resistance 
corresponding to a fixed edge is therefore zero, which means that 
any amount of current may flow without producing a potential 
difference. 


It is interesting to note that at a loaded edge or point the 
force against displacement characteristics cf the applied load 
may vary from a dead weight, where the force is independent of any 
motion of the end of the specimen, to a screw- type loading machine, 
where the applied force varies greatly with slight motions of the 
specimen. Between these two extremes would be a structure which is 
loaded by another, the stiffness of which is of a corresponding 
order of magnitude. - 
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Electrically, the preceding loading conditions correspond to 
current supplied "by a constant -current generator, a constant- 
voltage generator, and a generator with an internal resistance 
intermediate between infinity and zero. The present computer is 
built with constant-current generators and, therefore, obtains 
solutions which apply to the dead-weight type of loading since 
this type of loading agrees more closely with the loading of the 
actual structure . 

The differences between the computer solution, wherein the 
structure is broken up into a finilo 'lumber of sections, and the 
exact solution, wherein the number of sections is increased 
without limit, should be considered. 

In figure 7 a simple panel is shown loaded at one end of the 
flange. The axial force in the flange gradually "leaks off" as 
shear force into the sheet; the shear force in turn is transmitted 
as axial force into the stringer. The amount of force leaking off 
will be determined by the shear stiffness of the sheet and the 
axial stiffnesses of the flange and stringers . The smooth curve 
shown in figure 7 thus represents the stress in the flange. The 
curve is high at the loaded end and decreases smoothly as the 
distance from the end increases. 

In order to set up the equivalent electrical network, the 
structure must be broken into a finite number of sections so that 
a section of resistor network may be adjusted to correspond to 
each section of the structure. As the distance increases uniformly 
along the electrical network, the current in the resistors 
representing the flange remains constant in any one resistor and 
then abruptly decreases to the next lower value as a Junction 
point is reached. This condition is .graphically illustrated by 
the series of steps in figure 7* At the points where the step 
function crosses tho smooth curve there is, of course, no error 
due to the finite steps. At all other points a relative error 
will exist between the two curves) this error may generally be 
reduced in magnitude by adjusting the steps to lie half above 
and half below the smooth curve . 

In order to approach this condition, the stiffness of each 
section is considered to lie at its center, and the sections are 
shifted so that their centers come over the points Lq, and 

Lg in figure J. The first section must therefore be made of 
half length with the stiffness concentrated at the loaded end) 
the last section will also be of half length but with the stiffness 
concentrated at the point L^, or root. Since the end conditions 
fix the shear at the root as zero, this last half section does not 


MCA TN No. 1281 


7 


enter into the calculation. When the stiffnesses are adjusted in 
this manner, a smooth curve drawn through the midpoint of each step 
will closely approach the exact solution- 


NUMERICAL EXAMPLES OF SHEAR -LAG PROBLEMS 
Analysis of a Sheet -Stringer Panel 


As an illustrative example the panel shown in figure 8 was 
chosen for analysis- A single -stringer structure was used 
in order to compare results with the exact solution as obtained 
by the method given in appendix B of reference 2. 

As shown in figure 8, the numerical values of the symbols are 


b = 12 inches 
t = 0.024 inch 
Ap = 0.228 square inch 
A l = 0-500 square inch 
L = 100 inches 

Also 

E = 10,400 ksi 
G = 4000 ksi 

In order to agree with the configuration of the computer network, 
the structure was divided lengthwise into five equal parts- 
Substituting the foregoing values in equation (2) and choosing 
a = 70,000 results in flange resistances of 


= 70000 X 20 

F 10400 X 0-228 


= 590 ohms 


and in stringer resistances of 


p „ 70000 X 20 

L lo 4 oo x 0.500 


269 ohms 
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As explained previously, the length of the first shear hay 
will he one-half of the remaining hays, or 10 and 20 Inches, 
respectively. From equation 4, the shear resistances will 
therefore he for the first hay 


r-i 7.0000 x 12 s, 875 ohms 

4000 x 0.024 x 10 

and for each of the remaining hays 


r - 7C000 x 12 

2 4000 X 0.024 x 20 


=437 ohms 


The Input current I is chosen for convenience at 15, milliampere 
this value corresponds to the load of 1000 pounds. Then hy 
equation (1) 


(3 « _H_ = 0 .015 
1000 

These values of R, r, and p were set up on the computer 
network and the corresponding currents were recorded as shown in 
figure 9* Substitution of the current values in equations (l) 
and ( 3 ) then gave the desired stress values. For instance, the 
current in the first shear resistance r-j_ was 5 *90 milliamperes 
and hy equation ( 3 ) the shear stress was 


t = 119 . = 1640 psi 

0 .015 X 10 X 0 .024 

The results of the analysis are shown in the plots of figure 10. 
Normal stresses derived hy the computer method are shown as steps 
for reasons previously explained, with the computed exact-method 
curves superimposed. The exact-method curve of longitudinal stress 
does not cross the last step exactly in the middle because of the 
rapidly changing slope of the curve in this region. The plot for 
shear stress hy the computer method is shown as a set of "standpipes" 
inasmuch as the transverse current corresponding to the shear can he 
measured only at the transverse junction points. The maximum value 
of stress given hy this method does not agree very well with the 
maximum value of the exact method. 
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The results can he improved by adjusting the length of the hays 
that is, hy using shorter lengths in the regions of steeper slope. 

In order to note the effect of this adjustment, the computer was 
set up with hay lengths of 50, 20, 15, 10, and 5 inches; and a new 
analysis was made. The rosults are shown in figure 11. The 
agreement with the theoretical curve is now excellent, hoth for 
axial stresses and for shear stresses. 


Analysis of a Box Beam 

When the cover of a hox beam is to he analyzed, the resistances 
are deteimined and set up exactly as for the panel. The load, 
however, is now introduced into the flange hy shear forces in the 
weh. When the beam is loaded at the tip, the shear is uniformly 
distributed. In the network, however, it is, of course, necessary 
to introduce the current in discrete amounts at finite intervals. 

As for the previous case of the panel in which the shear resistances 
of the sheet were distributed, the current introduced at the free 
end of the hox beam is made one -half of the amount introduced at 
succeeding points so as to obtain a series of steps which straddle 
the smooth curve . 

As a numerical example the panel previously considered is 
used as the cover of a hox beam. (See fig. 12.) Beca\xse the 
beam is symmetrical about the center line, only one-half of the 
beam need be considered. The depth h is 3 inches, and a load 
of 250 pounds is assumed on each side at the tip. Also, the 
beam is divided into five equal-length bays of 20 inches. The 
running shear per unit length of the flange is 


P 250 

h" 3 


83.3 pounds per inch 


For a bay length of 20 inches the total shear force per bay is 
83.3 X 20 = 1666 pounds. As before, a convenient value of 
10 milliamperes is chosen for the current, and the constant p 
then is 


3 = — =2— = 0.0060 
1666 


The equivalent network is shown in figure 13, and the 
corresponding stress-distribution plots are shown in figure 14. 
Since the slope of the exaoc-method curves does not change very 


1C 


NACA TN No, 1281 


rapidly, the use of equal- length bays is shown to give very good 
agreement "between calculated and tost results. It should he noted 
that, with the structure divided into equal hays with a half step 
at the loaded, end, a half step Is left unaccounted for at the root 
end. This omission may he. disregarded without appreciable error. 

ANALOGY FOB BOLTED -JOINT PROBLEMS 


As shown In appendix A, the correspondence between the 
elements of a bolted Joint (fig. 15) and the computer network 
(fig. 16) is .given in the following table i 


Bolted Joint 

Electrical network 

Quantity 

Symbol 

Quantity 

Symbol 

f • * ; • . • #? * ^ 

Axial load 

P. 

Longitudinal current 

I 

.Axial flexibility 
of plate 

SK_ 

P 

Longitudinal resistance 
simulating plate 

R P 

Axial flexibility 
of strap. 

K 

s 

Longitudinal resistance 
. simulating, strap . 

*s 

Bolt, lead 

i. 

Transverse, current > 

i . 

Bolt constant 

. ,e . ... 

Transverse resistance 

r 

Axial displacement 

• > ’ *.• * • ! •“ 

u 

Potential 

V 


Vritten in equation form. with arbitrary constants a, and p, 
the following rolations hold: ' 


I = PP 
I'.m PB 

*p - 

S B - <*6 

r = oC 
7? flfiu 


( 6 ) 

( 7 ) 

( 8 ) 
(9) 
(1C) 

(ID 
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From a comparison of figures 15 and 16 , the resistor network 
will he noted to consist o^ a single row of longitudinal resistances 
E e which simulate the two straps of the "bolted Joint. A single 
set of transverse resistances r n which simulate the two shear 
areas of the bolt is also used. These apparent discrepancies are 
accounted for by the appearance of the factor 2 in equation (8). 

When bolted- Joint problems are to be solved, it should be 
noted that the load is transmitted through the bolts in discrete 

amounts acting at definite points, rather than through infinitesimal 
elements uniformly distributed along a sheet as in shear-lag 
problems. There is, therefore, no "finite-length increment" to 
cause errors when the network is used to simulate a bolted Joint. 


NUMERICAL EXAMPLE OF A BOLTED- JOINT PROBLEM 


The bolted Joint shown in figure 15 is considered with the 
following constants: for steel bolts, 


d = 0.25 inch 


E^ = 29,000 ksi 
and for plates of 24 S-T aluminun alloy, 


t = 0.3125 inch 

r 

t - O.I875 inch 

S 

p = 1.0 inch 


b a 2.0 inches 


E = 10,400 ksi 


The plate constants may be computed by the methods given in 
reference 3 * From equation ( 2 ) of reference 3 


K * _3L. 

8 bt a E 


2 x O.I875 x 10400 


1.0 


= 0.000256 


and 



= JP 


= 0.000308 

2 x 0.3125 x 10400 
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The holt constant is {riven "by equation (A19) of reference 3 



When t is replaced by t av , which is given by equation (A2l) 
of reference 3, as 


* = = g_X O.1873 + 0 » 3.125 ,0.344 inch 

aV 2 2 


the bolt constant becomes 


C 


8 


05 * 
P b 




8 

0.3125 x 29000 



+ 1 . 


0 .00229 


In setting up the resistance network, a convenient value for 
the constant a is 200,000; therefore, from equations (8), (9\ and 
(1C), respectively, 

Rp = 2C0000 x 0.000308 = 61.6 ohms 
= 200000 X 0.000256 =51.2 ohms 

8 

r =• 2C0000 X 0.00229 = 458.0 ohms 


These values were set up on the computer, and the transverse 
currents were measured and expressed as fractions of the total 
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applied current. For comparison, the corresponding "bolt loads 
were computed "by the methods of reference 3, and the results are 
tabulated as follows: 


Bolt 

Bolt load, B 
(fraction of total) 


Analytical method 

Electrical method 

1 

0.256 

0.253 

2 

.185 

.183 

3 

.159 

•159 

k 

.173 

.169 

5 

.223 

.228 


CONCLUSIONS 


The electrical -network computer using available variable resistors 
with sources of current stabilized by simple electronic regulating 
circuits is a practical device for rapidly solving problems in 
shear lag and bolted joints. Sufficient accuracy is obtained to 
warrant the use of the comouter in ordinary engineering stress 
analysis. The use of a more extended electrical network capable 
of simulating more complicated structures would appreciably 
reduce the time necessary for their analysis. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. February 14, I9L7 
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APPENDIX A 


DEVELOPMENT OF ANALOGY BETWEEN LOADS IN BOLTED 


JOINT AND CURRENT IN LADDER -TYPE NETWORK 


Figure l6 shows a simple electrical ladder-type network set 
up to simulate a "bolted joint. A current I is introduced into 
the network with part flowing down the line of resistances Rp 
and with part flowing through the transverse resistances r 
into the line of resistances R g . The resistances Rp are all 
equal, as are R 8 > "but the resistances r are not necessarily 
equal . 

Applying the principle that the voltage drop around any 
closed mesh must he eero results in 


in r n + I sn P 3 


in+1 r n+l " ^pn ^p = ^ 


(Al) 


or 


i-n+1 r n+l - in r n + ^an " ^pn ®p 


(A2) 


tut 


-XL. 



(A3) 


and 




(A4) 
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Substituting equations (A3) and (Ab) in equation (A2) therefore gives 


n 

%+! r n+l = "n ^ + (®g + Pp) 


Now 


n n-1 

> 1=1 + ~> 
— “ n 


m 


(A5) 

(a6) 


Substituting equation (A 6) in equation (A5) and rearranging terms 
results in 

F„ + R. F 

^ 4 - 

n 


4 _ x n 4 , A 's 

i n+ i = L + 


_P i 


n+1 


r n+l 


'n 


x i + 


n+l 


r n+l m 


(AT) 


From equation (l) of reference 3 the corresponding mechanical 
equation is 

C„ 2K_ + K 2 K £&_ + K a J3cl 

3 X ~ (A8) 


Bn*l - ®„ + - 1 - 7 " B „ - P * 


"n+l 


Cn+1 


c n+l 


11 M 

y~ b . 


"n+1 


m 


The following analogy is therefore obtained: 


Bolted joint 

Electrical network 

Quantity 

Symbol 

Quantity 

Symbol 

Axial load 

? 

Longitudinal current 

I 

Axial flexibility 
of plate 

pr 

P 

Longitudinal resistance 
simulating plate 

R P 

Axial flexibility 
cf strap 

K s 

Longitudinal resistance 
simulating strap 


Bolt load 

B 

Transverse current 

i 

Bolt constant 

C 

Transverse resistance 

r 

-Axial displacement 

u 

Potential 

V 
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Figure I- Simple stiffened panel. 
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Figure Z.~ Elementary bolted joint. 








Figure 3.- Front panel of electrical computer used in solution of 
shear -lag and bolted -joint problems. 
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Fig. 4 


“O Closed circuit jack 
jVW\T 1000- ohm adjustable resistor 
fc^l Constant current source (fig.5) 
Switch 





Figure 5- Wiring diagram of one constant 

current source. 
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Fig. 6 



Figure 6.- Rear view of panel of electrical computer used in solution 
of shear -lag and bolted -joint problems. 
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Fig. 7 
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Figure 7- Stiffened panel and 

equivalent network. 


Fig. 8 


NACA TN No. 1281 


P=IOOO 


t 


L= 100 


y//Y////////////////A'/ / 

— b = l2 


i— 1 = 0.024 
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Figure 8- Stiffened panel for 
numerical example. 
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Fig. 9 



Figure 9.- Equivalent network for panel 
showing magnitude and direction of 
current. RF=590ohms; RL=269ohms,- 
r, = 875 ohms; r 2 = 437 ohms. 


Normal stress, ksi Shear stress, ksi 


Fig. 10 
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Figure 10.- Stresses in panel of equal bays. 


Normal stress, ksi Shear stress, ksi 
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Fig. 11 



Figure 1 1.- Stresses in panel of unequal bays. 


Fig. 12 
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Figure 12.- Box beam for numerical example. 


NACA TN No. 1281 


Fig. 13 
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Figure 13.- Equivalent network for box beam 
showing magnitude and direction of current. 
R F =590ohms; R L =269ohms; r, =875 ohms; 
r 2 = 437 ohms. 
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Figure 14- Stresses in beam cover. 
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Figure 15- Typical bolted joint. 
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Fig. 16 
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Figure 16.- Equivalent network for 

bolted joint. 


